Phase diagram of quantized vortices in a trapped Bose-Einstein condensed gas 
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We investigate the thermodynamic stability of quantized 
vortices in a dilute Bose gas confined by a rotating harmonic 
trap at finite temperature. Interatomic forces play a crucial 
role in characterizing the resulting phase diagram, especially 
in the large N Thomas-Fermi regime. We show that the crit- 
ical temperature for the creation of stable vortices exhibits a 
maximum as a function of the frequency of the rotating trap 
and that the corresponding transition is associated with a dis- 
continuity in the number of atoms in the condensate. Possible 
strategies for approaching the vortical region are discussed. 

PACS numbers: 03.75.F, 05.30.J, 32.80 Pj, 67.40.Db 



The occurrence of quantized vortices in superfluids has 
been the object in the past of fundamental theoretical 
and experimental work jl]-|5| . After the experimental re- 
alization of Bose-Einstein condensation in trapped alkali 
gases ||, the structure of vortices in these novel many- 
body sytems has soon attracted the attention of theorists 
(for a recent review on various theoretical aspects 
of Bose-Einstein condensation in trapped gases, see (HI ) . 
Special attention has been devoted, in particular, to the 
calculation of the critical frequency needed to create a 
stable vortex in the rotating frame. This quantity is of 
crucial importance in view of the experimental possibil- 
ity of creating vortices by rotation of the confining trap. 
In systems interacting with repulsive forces the critical 
frequency turns out to be smaller than the oscillator fre- 
quency characterizing the confining harmonic potential 
and to decrease smoothly with the number of atoms in 
the condensate. The purpose of this work is to discuss the 
behavior of the critical frequency at finite temperature, 
pointing out the crucial role played by the interactions, 
especially in the large N, Thomas-Fermi regime, where 
the system behaves like a Bogoliubov gas. We will show 
that, due to the harmonic nature of the confinement, the 
phase diagram for the vortical configurations of these di- 
lute Bose gases exhibits new interesting features which 
are absent in dense superfluids. 

Let us consider a trap rotating with frequency u> along 
the z-axis. The equilibrium configuration of the sys- 
tem in the rotating frame is obtained by minimizing 
the free energy F = E — TS, using the Hamiltonian 
H(lo) = H — ldJz where H and J z are the Hamiltonian 
and the angular momentum (third component) of the 
many-body system, evaluated in the laboratory frame. 
A vortical configuration is thermodynamically stable if 



the trap rotates with frequency uj larger than the critical 
frequency 
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fixed by the difference between the free energy of the 
configurations with and without the vortex and by the 
corresponding difference of the angular momentum. The 
stability criteria for vortices, including the possible occur- 
rence of metastable configurations, have been the object 
of several recent theoretical studies jy] |l4| . In particular 
one has to distinguish between configurations where the 
system is stable against small deformations of the system 
(local stability) and configurations where the system is 
in full thermodynamic equilibrium (global stability) . For 
frequencies smaller than the critical frequency ([!]) the 
system, in the presence of a vortex, is not thermodynam- 
ically stable, but can be locally stable (metastability) 
pd] , ^2[ . For even smaller frequencies the vortex is unsta- 
ble also locally and exhibits elementary excitations with 
negative energies jl3|,[l5]], unless one considers toroidal 
geometries. In the following we will discuss only configu- 
rations in which the system is in full thermal equilibrium 
in the rotating frame. 

We will consider a gas confined by an axially symmetric 



harmonic potential of the form V e: 
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-y . Actually, in order to achieve thermal- 
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ization in the rotating frame the trap should contain an 
asymmetry in the xy-plane. One however expects that, 
unless the gas is too cold, even very small asymmetries 
are sufficient to ensure fast thermalization, so that one 
can safely use symmetric potentials for the calculation 
of the vortical configuration as well as of the relevant 
thermodynamic functions and, at the same time, assume 
thermal equilibrium in the rotating frame. 

In the absence of two-body interactions the critical 
frequency uj c is easily calculated and coincides with the 
radial oscillator frequency luj_. In interacting gases the 
value of u> c can be calculated, at zero temperature, by 
solving the Gross-Pitaevskii equation for the order pa- 
rameter H). In this case equation (||) takes the simpler 
form to c = (E v — Eq)/NTi where E v and Eq are the en- 
ergies of the ground state configurations with and with- 
out the vortex, respectively, and we have considered a 
vortex with one quantum of circulation: (J z ) v = Nh, 
(Jz)o — 0. Typical values of uj c correpond to a few tens 
of Hertz. Analogously, for each value of to, one can cal- 
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culate a critical value Nq such that, for Nq 



N§ the 
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vortex corresponds to a stable configuration in the frame 
rotating with frequency to. Actually the Gross-Pitaevskii 
equation exhibits a scaling behaviour in the interaction 
parameter Noa/aho, where a is the s-wave scattering 
length, which will be always assumed to be positive, and 
dho = y / njmujh is the oscillator length relative to the 
geometrical average ujho = (^z^j) 1 ^ 3 = A 1 / 3 ^ where 
A = w 2 /wi is the deformation parameter of the trap. As 
a consequence of scaling one can always write the critical 
number An in the form 



N§(u>)- 
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where f\ is a function that can be calculated by solv- 
ing the Gross-Pitaevskii equation for each value of A. 
When the dimensionless parameter N^a/aho is much 
larger than 1 (Thomas-Fermi limit) the solution of the 
Gross-Pitaevskii equation gives the following result |l6[ ] 
(see also for the critical frequency 
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where R±_ = A 1,/3 a^ (IbN^a/ aho) 1 ^ is the Thomas- 
Fermi radius of the condensate in the xy-pleme, while 
£ = a ho (15A a/a/ lo ) -1 ^ 5 is the healing length, which 
fixes the size of the vortex core. The Thomas-Fermi 
result ([5]) allows one to determine in a simple way the 
function fx which turns out to scale as A -5 / 6 . In the op- 
posite limit, when the parameter Noa/aho is very small, 
one can use a perturbative treatment of the interaction 
yielding the result uj c /u±_ = 1 — A 1 / 3 (87r) _1 / 2 Aoa/a? lo . 
The quasi ideal gas limit is not however relevant for the 
available experimental configurations since N^a/aho is al- 
ways much larger than 1. In fig. 1 we compare the results 
for the function A 5 / 6 /a, calculated by solving the Gross- 
Pitaevskii equation with the prediction of the Thomas- 
Fermi approximation for two different values of the de- 
formation parameter (A = \/8 and A = 1). One sees that 
the Thomas- Fermi limit is rather accurate for u> < 0.5wj_. 

From an experimental point of view the creation of a 
vortex by adiabatic increase of the rotational frequency 
of the trap up to u> c might not be the safest procedure 
at very low temperature. In fact one expects that in this 
case the system will exhibit a high barrier for the nu- 
cleation of vortices as happens in superfluid helium ||. 
Though the nucleation process in dilute gases might be- 
have differently from the one of a dense supefluid, one 
nevertheless expects that a safer procedure to generate 
vortices is first to rotate the gas at higher temperature 
(in case, also above the critical temperature for Bose- 
Einstein condensation) and then to cool it via evapora- 
tion. This strategy has motivated the present investiga- 
tion of the phase diagram of vortices at finite tempera- 
ture. 



Let us first discuss the behaviour of the critical tem- 
perature T c as a function of the rotational frequency. 
We will suppose the system to be in thermal equilib- 
rium in the frame rotating with frequency u. A first 
estimate of the critical temperature is obtained by tak- 
ing the ideal gas model. In the absence of rotations 
and for large values of N this model predicts the result 
T c ° = 0.94(^ /lo /fc s )A 1 /3. 

The effect of the rotation is to change the shape of the 
distribution function of the thermal atoms so that the 
number of atoms out of the condensate takes the form 
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Result (^) holds for T < T c where the chemical potential 
of the non interacting model vanishes in the thermody- 
namic limit. Notice that, unless a vortex is formed, the 
rotation has instead no effect on the condensate, whose 
wave function has zero angular momentum since the con- 
fining trap is axially symmetric. By using the transfor- 
mation p —>■ p — mu) x r the integral (Q) takes the same 
form as in the absence of rotation, with an effective radial 
frequency luj_(lu) = {uj 2 ^ — oj 2 ) 1 / 2 . This reflects the role 
of the centrifugal force which is consequently responsible 
for a lowering of the critical temperature T c according to 



T c (io) = T C U 1 - -y 
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where T~ is the critical temperature in the absence of 
rotation. In the same model the temperature dependence 
of the condensate fraction is given by: 
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The lowering of T c is dramatic near ui± where the sys- 
tem becomes unstable. In this case the centrifugal force 
exactly cancels the confining effect of the harmonic po- 
tential. For confining potentials growing faster than r 2 at 
large distances the critical temperature would instead dif- 
fer from zero for any value of the rotating frequency. The 
behaviour exhibited by dilute gases in harmonic traps 
should be compared with the one of dense liquids, like 
superfluid helium, where the critical temperature is prac- 
tically unaffected by the rotation. 

The next step is to assume that the T = relation- 
ship Nq(u>) of eq.(||) between the frequency uj and the 
minimum number of atoms in the condensate required to 
have a stable vortex, is not changed at finite temperature. 
This relationship allows one to calculate, through eq.(||), 
the critical temperature T v , below which the vortex cor- 
responds to a stable configuration in the frame rotating 
with frequency to. We find the useful result 
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where we have explicitly used eq. (||) for the critical tem- 
perature. Equation (|7|) points out the crucial role played 
by two-body interactions. In particular, by increasing 
the value of the dimensionless parameter Na/a,ho, the 
curve T v is pushed to the left, making the region of sta- 
bility for vortices wider and wider. In fig. 2 we show the 
critical curves T c (uj) (dashed line) and T v {tu) (full line). 
The function fx entering eq.(|^) was obtained by solving 
the Gross-Pitaevskii equation. For temperatures below 
T c (lo) the gas exhibits Bose- Einstein condensation. How- 
ever only for temperatures below T v (lo) will the vortex 
correspond to the equilibrium configuration of the sys- 
tem. 

In deriving the main result (Q) we have made sev- 
eral approximations that it is worth discussing explic- 
itly. These concern the use of the ideal gas predictions 
(||) and (JsJ) for the critical temperature and the thermal 
depletion of the condensate, as well as the use, at finite 
temperature, of the T = relationship between the ro- 
tational frequency and the critical value of Nq, given by 
the solution of the Gross-Pitaevskii equation. 

Neglecting interaction effects in the evaluation of the 
critical temperature is a very good approximation since 
near T c the gas is extremely dilute. Indeed self-consistent 
mean field (110 as well as ab initio Jl9| calculations 
show that T c is lowered by interactions only by a few per- 
cent. The effects on the thermal depletion of the conden- 
sate are instead expected to be more important because 
of the sizable renormalization of the chemical potential 
due to the presence of the condensate. For the same 
reason one expects that the relationship (|J) between the 
frequency and the critical value of Aq should exhibit a 
temperature dependence. 

In order to understand the importance of the above ef- 
fects let us calculate the chemical potential /i = dF/dN 
of the gas with and without the vortex. Also at fi- 
nite temperature the chemical potential can be calcu- 
lated, with good approximation, by solving the Gross- 
Pitaevskii equation for the condensate with a fixed value 
of Nq. In the absence of vortices and for large val- 
ues of N^a/aho one can use the Thomas- Fermi formula 
Ho(N ) = ^huiho(^5N a/aho) 2 ^ 5 ■ If a vortex is present 
the chemical potential will contain an extra term <5/i = 
!?kj c with uj c given by eq.(^). Notice that both /io and 
Sfi are fixed by the value of No. However, for a given 
value of N and T, the value of Ao differs depending on 
whether the vortex is present or not. This effect is absent 
if one uses the ideal gas prediction (^). By treating the 
interaction effect as a small perturbation one finds the 



result 1 10 1 



which generalizes the ideal gas prediction (||) ■ The above 
equation shows that, because of interactions, the num- 
ber of atoms in the condensate with and without the 
vortex differ, to the lowest order in the interaction, by 
the amount 
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The discontinuity 5 No in the number of atoms in the con- 
densate is not very large (at the maximum of the curve 
T v of fig. 2 one finds 5N ~ 200), but should be taken 
into account for a safe estimate of the critical frequency. 
In fact the relevant change in the chemical potential is 
given by fx v (N,T) - fi (N,T) = (dfi /dN )SN +Sfi. By 
using (g) and @) one finds that the leading corrections, 
due to the interaction term in (|J), cancel out. This sug- 
gests that, at finite temperature, the critical frequency 
can be safely calculated using the T = solution of the 
Gross-Pitaevskii equation and the ideal gas expression 
(H) for Ao In conclusion the main result (Q) for the crit- 
ical temperature for vortices is expected to represent a 
rather good estimate also when T v is comparable to T c . 
Notice that the value of Ao on the line T v can be never- 
theless significantly smaller than the ideal gas prediction 
(H) because of the interaction term in (ph. This effect 
should lower the height of the barrier at the transition 
and hence favour the nucleation of the vortex. 

A possible way to approach to the vortical region 
is through an adiabatic transformation. Suppose that, 
starting from a system initially at rest in the labora- 
tory, we switch on the rotation of the trap in an adi- 
abatic way. Since the entropy of the gas is given, 
in first approximation, by the ideal gas value S — 
4fc s A(C(4)/C(3))(T/T c (w)) 3 the adiabatic increase of the 
frequency has the effect of lowering the temperature. If 
the initial temperature is below T c the transformation 
will eventually bring the system into the vortical region. 
Notice that this scenario does not require that the initial 
temperature be much smaller than T c . This is impor- 
tant for at least two reasons: first the presence of a large 
fraction of thermal atoms is crucial to favour the ther- 
malization of the gas in the rotating frame; second, the 
nucleation of the vortex may become easier if the number 
of atoms in the condensate is not too large. 

Once the vortex is created a major problem is its ex- 
perimental detection. Recently, various methods have 
been proposed to detect vortices in these trapped gases. 
They include the imaging of the vortex core following the 
expansion of the atomic cloud pQ| , the occurrence of fre- 
quency shifts in the collective oscillations and the 
occurrence of dislocations in the interference patterns be- 
tween condensates containing vortical configurations p2~| . 

In the most interesting vortical region the system is ex- 
pected to exhibit further interesting features. In fact in- 
creasing the rotational frequency will favour the creation 
of more complex vortical configurations, associated with 



3 



the occurrence of 2 or more vortices |2^;|24j . In principle 
vortices with higher quanta of circulation are also pos- 
sible. However these configurations are likely unstable 
Hflij. A final interesting question concerns the stabil- 
ity of elementary excitations in the rotating gas. In ]2j| 
it has been shown that, for large frequencies, the surface 
excitations of the condensate become unstable. This hap- 
pens for frequencies larger, but not tremendously larger, 
than the critical frequency needed to create a vortex. 
For example with TV = 10 4 atoms in a symmetric trap 
at T — it is found |^5| that for the critical frequency 
for generating a surface instability is lu/loj_ ~ 0.5, while 
the one for the vortex instability is ~ 0.35. The influence 
of the vortex on the stability of such excitations jl2| as 
well as the possible occurrence of new scenarios associ- 
ated with the "condensation" of surface excitations are 
challenging questions for future investigation. 
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Figure Captions. 

Fig.l. Function A 5 / 6 f\ calculated for two different val- 
ues of A (open circles: A = 1; solid circles: A = The 
prediction obtained from eq.(|^) is also reported (dashed 
line). Frequencies are in units of luj_. 

Fig. 2. Phase diagram for vortices in a harmonically 
trapped Bose gas (N — 10 4 , a/a/ lo = 7.36 x 10~ 3 and 
A = 1). Frequencies are in units of u±. 
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